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Abstract 

We study the average density of resonances, (p(x,y)}, in a semi-infinite dis- 
ordered chain coupled to a perfect lead. The function (p(x,y)) is defined in the 
complex energy plane and the distance y from the real axes determines the reso- 
nance width. We concentrate on strong disorder and derive the asymptotic behavior 
of (p(x, y)) in the limit of small y. 



1 Introduction 

Open quantum systems often exhibit the phenomenon of resonances. Resonances corre- 
spond to quasi-stationary states which have a long life-time but eventually decay into the 
continuum. (A particle, initially within the system, escapes to infinity.) They are charac- 
terized by complex energies, E a = E a — |r a , which correspond to poles of the S-matrix 
on the unphysical sheet of the complex energy plane [1, 2]. There are many examples of 
resonances in atomic and nuclear physics. More recently, there has been much interest in 
resonant phenomena in the field of chaotic and disordered systems (for recent reviews see 
[3,4]). 

There is considerable amount of work concerning the distribution P(T) of resonance 
widths in one-dimensional disordered chains [5, 6, 7, 8, 9]. Numerical studies presented 
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in that work demonstrate that, in a broad range of T, P(T) oc T" 7 , with the exponent 7 
being close to 1. (This behavior is not restricted to disordered chains, but pertains also to 
two- and three-dimensional systems with localized states [7, 9]). An intuitive argument, 
which assumes a uniform distribution for the localization centers of exponentially localized 
states, leads to a (1/T) -behavior [7, 8, 9]. The analytical calculation in [8], performed for a 
one-dimensional continuous (white-noise) potential, exhibits this behavior for sufficiently 
small T, followed by a sharp cutoff at still smaller T, due to the finite size of the sample. It 
is also shown in [8] that in a broad range of T, -P(r) is well fitted by a function T" 1 - 25 . In 
the present Letter we develop an analytical approach for a discrete, tight-binding random 
chain. We treat the problem in the limit of strong disorder and derive the asymptotically 
exact (1/T) -behavior for a semi-infinite system. 

2 The model and its Effective Hamiltonian 

We consider a semi-infinite disordered chain coupled to a (semi-infinite) perfect lead. 
n = 1,2,...., designate sites along the chain. Each site of the chain is assigned a site 
energy, e n . Different e ra 's (n = 1,2,....) are independent random variables chosen from 

some symmetric distribution q(e). Sites of the lead are labelled by n = 0, —1, —2, All 

sites of the lead are assigned e n = 0. The lead simulates the free space outside the chain. 
All nearest neighbor sites of the full system (chain + lead) are coupled to each other by 
a hopping amplitude t, so that a particle, initially located somewhere within the chain, 
will eventually escape into the lead. 

The most direct approach to the problem of resonances amounts to solving the sta- 
tionary Schrodinger equation, for the entire system, with the boundary condition of an 
outgoing wave only. This condition, which makes the problem non- Hermit ian, describes a 
particle ejected from the system. The Schrodinger equation with such boundary condition 
admits complex eigenvalues E a , which correspond to the resonances [1, 2]. This kind of 
approach, which leads in a natural way to a non-Hermitian effective Hamiltonian, has 
been used for a long time in scattering theory, including scattering in chaotic and disor- 
dered systems ([4, 10], and references therein). For our system the approach amounts to 
solving the infinite set of coupled equations 

-tlp n+1 - tljln-! + e n ip n = Eip n (-00 < n < 00) (1) 

We recall that e n = for n < 1 (the lead) and it is random for n > 1 (the chain). Eqs.(l) 
are to be solved subjected to the boundary condition ip n oc exp — ikn , corresponding to 
an outgoing wave in the lead. The complex wave vector k determines E according to 
E = —2tcosk. It is straightforward to eliminate from Eqs.(l) all ^„'s with n < 1 , thus 
reducing the problem to a system of equations for the amplitudes ip n on the sites of the 
disordered chain alone (n = 1, 2, ....). 

-t^ n+1 - tyn-x + l n i) n = Ei/; n (n = 1, 2, . . .) (2) 
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with the condition ipo = 0. Here e n = e n for n — 2, 3, but not for n — 1. This end site 
is assigned a complex energy l\ — e\ — t exp ik which describes coupling to the outside 
world. Thus, the resonances are given by the complex eigenvalues of the non-Hermitian 
effective Hamiltonian defined in (2). Note that Eq. (2) does not constitute a standard 
eigenvalue problem because 1\ contains k, which is related to E. The problem is often 
simplified by fixing k at some real value, consistent with the value of energy near which 
one is looking for resonances. Such simplification corresponds to replacing the "exact 
resonances" by "parametric" ones [5]. There are indications that for sufficiently narrow 
resonances the "parametric" and the "exact" distributions are close to one another. In 
this Letter we restrict ourselves to "parametric" resonances. For instance, close to the 
middle of the energy band we set k = ir/2, thus arriving at the effective Hamiltonian 

(Heff)nm = ^n&nm ~ t nm (n = 1, 2, . . .) . (3) 

where t nm = t for nearest neighbors (and zero otherwise), and e\ — e\ — it. Let us repeat 
that all site energies, except for ei, are real. The imaginary part —it of 1\ accounts for the 
coupling of the chain to the lead, via the hopping amplitude t connecting site 1 to site 0. 
It is convenient to slightly generalize the model by assigning to this particular amplitude 
a value t', which can differ from all the other hopping amplitudes t. This allows us to tune 
the coupling from t' — (closed chain) to t' — t (fully coupled chain). In what follows we 
set t — 1 and denote the dimensionless coupling strength v = t'/t. 



3 The Average Density of Resonances 

Resonances correspond to complex eigenvalues of H e ff, i.e., to the poles of the resolvent 

6 = V^H- (4) 

in the lower half of the complex energy plane. We denote these poles by z a = x a + iy a . 
(These are just the complex energies E a = E a — |r a , of Section 1, in units of t which we 
set to 1.) The average density of these poles is defined as 

(p(x,y)) = (%25(x-x a )8(y-y a )), (5) 

a 

where (• • •) denotes disorder averaging, i.e., averaging over all realizations of the set {e n } 
of the random site energies. 

In order to appreciate the difference between (p(x, y)) and the probability distribution, 
P(x, y), of resonance width y (for some fixed x) let us consider for a moment a finite chain, 
of iV sites. In this case P(x, y), by definition, is equal to {p(x, y)) divided by N . When N 
increases, P(x,y) "runs away" towards smaller and smaller y's, and in the iV — > oo limit 
it approaches S(y). (Indeed, for a semi-infinite chain an eigenstate will be localized, with 
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probability 1, at infinite distance from the open end and, thus, will be ignorant about the 
coupling to the external world.) On the other hand, (p(x,y)) does have a well defined 
N — > oo limit, for any fixed (non-zero) y. To clarify this assertion, we begin with a closed 
semi- infinite chain (v = 0). In this case all states are localized (y a = 0)- some close to the 
end (site l)and some further away. When the chain is opened (v — 1), localized states turn 
into resonances. The point is that for fixed y (and x) the main contribution to (p(x,y)) 
will come from states that (for v — 0) were localized around some optimal distance d(x, y) 
from the open end. Therefore distant pieces of the chain do not contribute to (p(x,y)) 
and a well defined N — > oo limit exists. The N necessary for achieving this limit will, of 
course, depend on y (and x) but the limit will eventually be achieved for any y, however 
small (different from zero). Thus, although for any finite N, P(x,y) and (p(x,y)) differ 
only by the normalization factor N, it is (p{x,y)) that has a meaningful N — > oo limit. 
(p(x,y)) may be expressed in terms of the resolvent G as [11, 12] 

(p(x, y)) = ^(d x + id y ) (TrG(x, y)) . (6) 

One can interpret the z a 's as the positions of unit electric point charges in the plane, 
which give rise to an electric field E = E x k + E y y. In this picture, one has simply 
TrG(:r, y) = E x — iE y , and (6) is then interpreted as the Poisson equation for the averaged 
electric field and charge density. Equation (6) holds also for a closed chain, if G(x,y) is 
replaced by G(x,y), where the untilded G is the resolvent of the Hermitian Hamiltonian 
of the closed chain. In this case all the charges must lie on the real axis. Therefore, for 
any y different from zero, one can rewrite (6) as 

(p(x,y)) = ±(d x + id y )(Tr(G-G)). (7) 

The advantage of this representation is that, although both TrG and TrG diverge in the 
N — > oo limit, Tr(G — G) remains finite. (Of course, it is still perfectly all right to use (6), 
with the proviso that the derivative with respect to y is taken before taking the N — > oo 
limit.) 

The electrostatic analogy just mentioned makes it clear that (TrG(x, y)) cannot be a 
complex analytic function of z = x + iy, but must depend separately on x and y. This 
lack of analyticity makes the computation of (p(x,y)) a more difficult task than that of 
the density of states for a closed (Hermitian) system. In the latter case there is a well 
known representation of the resolvent 

where Tl n (z) is the self-energy at site n. Let us stress that (8) applies to a specific 
realization of random site energies {e n }, i.e., Y, n (z) is a random quantity which depends 
on the set {e n }. The statistical treatment of Y, n (z) forms the basis of many studies of the 
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localization problem [13], starting with the original work of Anderson [14]. The notion of 
self-energy can be generalized to our non- Hermit ian problem, and the analog of (8) is 



Gnn(z) = ~-T- • (9) 

z — e n — l^ n \z) 

This equation is essentially a definition of T, n (z). It is convenient to assign the imaginary 
term —iv to rather than to the "bare" energy of the first site, so that all site 

energies, {e n }, in (9) are real and are the same as for the corresponding closed chain. The 
term — iv will appear as a boundary condition for the self energy and will play a crucial 
role, as demonstrated below. 

In terms of the locator expansion [13, 14] S„(z) can be represented as a sum of paths 
which start at site n, visit other sites and then return (only once!) to the starting point n. 
Since the paths consist of steps connecting nearest neighbors (to the left or to the right), 
it is clear that £ n (z) can be decomposed into two pieces, "left" and "right": 

t n {z) = L n (z)+Rn(z) (10) 

The "left" self-energy, L n (z), depends only on the energies of sites to the left of site n, 
i.e., on e/s with 1 < j < n. Similarly, R n depends only on e/s with j > n and, thus, 
it is ignorant about the fact that the chain is coupled to the outside world, via site 1 
(that is why the R n are not tilded). The reason for decomposing S n (z) into L n (z) and R n 
is that these quantities obey simple recursion relations, which can be iterated to obtain 
their probability distributions: 

L ] {z) = [z-e ] . 1 -L 3 . l {z)]- 1 (11) 

and 

R J (z) = [z-e J+1 -R J+1 (z)]- 1 (12) 

Relation (11) has to be iterated starting with j = 2, with the "initial condition" Li = —iv. 
Relation (12), for the semi-infinite chain considered in this paper, leads to a stationary, 
n-independent distribution for the variable R n . We will not need this distribution in the 
forthcoming calculation, restricted to the case of strong disorder. 



4 Strong Disorder 

Strong disorder means that a typical value of e n is much larger than t — 1, i.e., the 
distribution of site energies, g(e), is very broad. In this case the above recursion relations 
simplify considerably. The real part of the self-energy can be neglected in comparison 
with e„, so that at any site n the real part of the resolvent is ReG nn (z) (x — e n ) _1 , 
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which is the same as for the closed chain. Thus, the real parts of the two resolvents in 
(7) cancel out and one obtains 

(p(x,y)) = -±-d y (TrIm(G-G)). (13) 

Furthermore, ImG nn (see (9)) is dominated by ImL n (z) and can be approximated as 

ImG nn = —— — J " V ( 2 ~ 7rS(x - e n )sign(I n - y) (14) 

where, in order to avoid cluttering the notations, we denote ImL n (z) = /„. Since e n and 
I n are independent random variables, the averaging of (14) is done in two steps: first over 
e n , which gives the strong disorder limit of the density of states, q(x), and then over I n , 
with the yet unknown distribution P n (I). The latter is governed by the recursion relation 
(11), which in the strong disorder limit simplifies to 

'> = (15) 

This equation, as well as the previous one, applies for n > 2. Site n = 1 provides the 
initial condition I\ = —v. The recursion relation (15) "propagates" the information from 
the open end of the chain (site 1) to distant sites. Site 1 is special and should be excluded 
from the trace in (13). This is because I\ has a fixed value, different from y. Using 
the electrostatic analogy, one can say that a charge, located away from the point (x, y) 
cannot produce a singularity at this point . Actually, (15) is valid only up to a point 
when / becomes of order y ( recall that we are interested in y fixed but arbitrarily small, 
i.e., in the "small-?/" asymptotic behavior of (p(x,y))). When such a small value of I is 
reached, the initial condition, i.e., the value of v, is forgotten, and there is no contribution 
to (p(x,y))) (more precisely, both ImG and ImG become very small and, moreover, 
cancel each other). In summary, averaging (14) with the distribution P n (I) and taking 
the derivative with respect to y, one obtains from (13) (with the G-term neglected): 

oo 

(p(x,y)) = q(x)J2Pn(y)- (16) 

n=2 

Eq. (16), supplemented by the recursion relation (15), enables us to obtain the asymp- 
totic (1/y) -behavior, as we now explain. Note that the random variables /„ are negative, 
so that resonances are located in the lower half-plane of the complex variable z. It is more 
convenient to work with the absolute value, \I n \, and to define y as positive, i.e., (p(x,y)) 
is the average density of resonances with width y (at energy x). Defining, instead of \I n \, 
a new variable, t n = ln\I n \, and designating its distribution as W n (t), we have, instead of 
(16), 

1 oo 

{p(0,y)) = q(0)-Y,W n (t), t = lny, (17) 

V n=2 
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where we set x = (middle of the band). It follows from (15) and the initial condition 
\h\ — v — 1 that 

n—1 

* n = -2X>k;l> (18) 

i=i 

so that 

c+ioo c+ioo 

W B (f) =< / ^e p(t - tn) >= J ^V^p), (19) 

c—ioo c—ioo 

where c is some real number and 

v(p) = J deq(e)e 2pln ^. (20) 
Summing the geometric series in (17) we obtain: 

1 C+ /°° dp ^ pt _ is(p) 



= J tfT^y t = ln «- (21) 



C—IOO 



This equation makes sense if \v{p)\ < 1. A sufficient condition for this inequality is that c 
be negative and small, so that we have to choose the contour somewhat to the left of the 
imaginary axis in the complex p plane. In order to prove the (l/y)-behavior, in the small 
y limit, we have to verify that the integral in (21) approaches a constant (non-zero) value 
for t — > — oo. It is intuitively quite obvious that in this limit, and for |c| <C 1, the integral 
will be dominated by the vicinity of the point p — 0, near which v(p) pa 1 + 2p < In \ e \>. 
For strong disorder, < In \ e \> is positive and equal to the inverse localization length l/£ 
(in the middle of the band). Replacing in (21) (1 — v{p)) by (—2p/£), we obtain: 

c—ioo 

Note that, since t is negative, the integration contour should be closed in the right half- 
plane, where Rep is positive, so that an extra minus sign is acquired. The formal proof 
of equivalence between (21) and (22), in the t — > — oo limit, is achieved by studying the 
difference between the two expressions. Let's denote by J and J' the integrals in (21) and 
(22), respectively. In both integrals c cannot be set equal to zero, because the integrand 
would diverge at p — 0. However, this divergence is cancelled in J — J' , so that for this 
quantity one can integrate directly along the imaginary p-axis. Moreover, since v(p), for 
p purely imaginary, is just the characteristic function of the distribution for the variable 
2ln\e\ : it must decrease faster than l/\p\ for large values of p. Therefore one can use the 
Riemann-Lebesgue lemma to prove that J — J' approaches zero in the t — > — oo limit. 
This completes the proof of the main result, (22), of our paper. Equation (22), gives the 
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asymptotically exact expression for the average density of resonances in a semi-infinite, 
strongly disordered chain. 
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